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摘要 :【 目的 】 荔 枝 蒂 星 吕 Conopomorpha sinensis Bradley 是 专 一 性 危害 我 国 华南 地 区 荔枝 和 龙眼 的 

重要 害虫 ,隐蔽 性 强 , 防 治 困 难 , 基 因 组 信息 缺乏 。 本 研究 的 目的 是 获得 荔枝 带 蛙 虫 的 基因 数据 , 导 

求 有 效 控 制 害 虫 的 分 子 靶 标 。 不 用 新 一 代 高 通 量 测序 技术 Ilumina HiSeq™ 4000 aF 25 4i e 

kk d SEAT EK ZU Jo M48 S00 2E, LR] EUR PEDE 68 996 条 unigenes。 进 一 步 利 用 

七 大 公共 数据 库 进 行 同 源 比 对 ,注释 了 22348 unigenes。 注 释 到 Nr 数据 库 的 unigenes 数量 最 多 ， 
达 27.0196 ,其 中 Nr 注释 的 荔枝 蒂 星 贝 unigenes P 5 4 X 9& Plutella xylostella unigenes 同 源 ， ids 

高 , 达 34.196 , 3$ unigenes 5 GO 数据 库 比 对 发 现 ,15 585 条 unigenes 根据 其 功能 大 致 可 分 为 3 类 

47 X, KEGG pathways 分 析 表 明 ,7 272 条 unigenes 定位 为 267 个 代谢 通路 。 基 因 注 释 进 一 步 得 

选 鉴定 获得 100 个 荔枝 带 星 虫 嗅觉 相关 基因 ;与 鲜 翅 目 相关 气味 结合 蛋白 基因 联合 分 析 发 现 ,与 萝 

枝 蒂 昱 忠 气 味 结合 蛋白 基因 直系 同 源 的 有 18 组 ,部 分 基因 形成 独立 一 比 。【 结论 ] 本 研究 首次 获 

得 了 荔枝 还 星 忠 的 转录 组 数据 ,研究 结果 为 生物 控制 iie aep. 了 重要 的 基础 数据 和 候选 分 

子 靶 标 。 荔 枝 东 星 贝 独 有 的 气味 结合 蛋白 基因 可 能 与 其 生境 中 特有 的 化 学 物质 相关 。 

关键 词 : DAER; 转录 组 ; 高 通 量 测序 ; 、 嗅觉 基因 
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Analysis of the  transcriptome and  olfaction-related genes of 
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Conopomorpha sinensis Bradley ( Lepidoptera: Gracilariidae) 

MENG Xiang"", HU Jun-Jie', LIU Hui', OUYANG Ge-Cheng', GUO Ming-Fang' ( 1. Guangdong 
Public Laboratory of Wild Animal Conservation and Utilization, Guangdong Key Laboratory of Integrated 
Pest Management in Agriculture, Guangdong Institute of Applied Biological Resources, Guangzhou 
510260, China; 2. School of Life Sciences, Guangzhou University, Guangzhou 510006, China) 
Abstract: [ Aim] Conopomorpha sinensis Bradley is a major and specific fruit borer pest of litchi and 
longan in Southern China. lt is difficult to be controlled because of its better concealment, and there is a 
lack of its genomic information. The objective of this study is to obtain the genetic data of C. sinensis and 
to seek effective molecular targets for pest control. [Methods] The transcriptome of C. sinensis was 
sequenced using the Illumina HiSeq"" 4000 sequencing platform and bioinformatics analysis. [ Results] 
The clean reads were then de novo assembled into 68 996 unigenes. Through a similarity search against 
seven public databases, 22 348 unigenes were annotated. Most unigenes were annotated to Nr database 
(27.0196) , and the unigenes of C. sinensis have the highest homology to those of Plutella xylostella in 
species distribution of unigenes in Nr database (34.196 ). According to GO database, 15 585 unigenes 
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were broadly divided into 3 categories and 47 sub-types. By searching KEGG database, 7 272 unigenes 


were located into 267 metabolic pathways. By further screening and identification, 100 olfaction-related 


genes of C. sinensis were obtained. Analysis of odorant binding protein genes from C. sinensis ( CSOBPs) 


and the closely related species of Lepidoptera revealed that these genes form 18 OBP ortholog groups, and 


some form a unique cluster. [Conclusion] This study acquired the transcriptome data of C. sinensis for 


the first time. The results will provide valuable basic data and candidate molecular targets for biological 


control of C. sinensis. The unique CsOBPs might be associated with habitat-specific chemical components 


of C. sinensis. 


Key words: Conopomorpha sinensis; transcriptome; high-throughput sequencing technology; gene 


annotation; olfactory gene 





35 BANJE HR Zee FR Es e p 0b. XC C5, R0 AR: to e o 
ART , 3 js E HR. Conopomorpha. sinensis Bradley ( f$ 
却 目 : 细 蛾 科 ) die — HE ft E xx BCRUJE BR B5) i R 
Hf (Thanh et al., 2006; 王 少 山 等 ,2008 )。 其 幼虫 
钻 蛙 为 害 嫩 梢 ` 嫩 叶 、 幼 果 和 成 果 , 运 用 化 学 防治 很 
难 奏效 ,并 且 致 使 果品 的 农药 残留 超标 ( 易 干 军 等 ， 
2002; 陈 炳 旭 等 , 2010) ,严重 影响 功 校 和 龙眼 的 食 
用 和 商品 价值 ,造成 重大 的 经 济 损失 。 随 着 我 国 逐 
渐 将 食品 安全 提 上 日 程 ,减少 农药 化 肥 的 使 用 ,做 到 
安全 无 公害 甚至 是 绿色 生产 ,使 消费 者 吃 到 放心 的 
果品 已 成 为 生产 优质 水 果 的 一 个 重要 内 容 。 其 中 ， 
33 Bos b B6 BEGUN GNO SS BUE" ESRTE 
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RNA 调控 等 问题 提供 了 很 好 的 发 展 前 景 ,并 迅速 成 
为 研究 鳞 翅 目 非 模式 昆虫 的 先进 技术 (杨帆 等 ， 
2014) 。 目 前 ,昆虫 嗅觉 基因 的 高 度 变 异性 使 得 同 
源 克隆 、 探 针 调 取 等 方法 无 法 获取 新 的 基因 序列 , 极 
大 地 限制 了 嗅觉 系统 研究 。 高 通 量 测序 技术 的 应 用 
对 于 推动 昆虫 嗅觉 的 研究 具有 重要 的 意义 。 

本 研究 以 TIumina Hiseq'" 4000 高 通 量 测序 技 
术 对 荔枝 带 蛙 虫 进行 转录 组 测序 研究 ,由 在 获得 更 
多 鳞 翅 目 非 模式 生物 功 枝 蒂 星 虫 的 转录 本 和 更 为 全 
面 的 转录 组 信息 。 通 过 发 展 高 灵敏 度 、 高 通 量 的 嗅 
觉 基因 识别 方法 ,从 转录 组 数据 中 发 掘 荔枝 蒂 星 虫 
生长 发 育 和 行为 调控 过 程 中 的 重要 基因 信息 ,并 与 
























































解决 的 重要 问题 。 

生物 防治 是 应 用 较为 广泛 的 无 公害 防治 技术 之 
一 ( 任 顺 祥和 陈 学 新 , 2012). HRT, zx SC RE HUE 
物 防 控 技术 研究 多 集中 在 天 敌 防治 植物 保护 剂 防 
治 和 攻 校 带 星 虫 挥发 性 的 化 学 生态 调控 物质 初 得 
( 王 少 山 等 , 2010; 黄 立 兰 等 , 2010) ,或 是 利用 初级 
合成 的 信息 素 干 扰 和 引诱 荔枝 蒂 星 虫 ( 彭 海 辉 等 ， 
2006, 2007) 。 然 而 ,生物 防治 对 功 校 带 旺 虫 作用 的 
分 子 台 标 尚 不 清楚 。 荔 枝 带 蝗虫 专 一 性 危害 荔枝 和 
龙眼 的 生物 学 特性 是 昆虫 对 寄主 在 长 期 进化 过 程 中 
相互 协同 进化 的 结果 (Karban and Baldwin, 1997; 
钦 俊 德 和 王 琛 柱 ，2001; Will and van Bel, 2006; 
Walling, 2008) 。 其 中 ,高 度 灵 人 敏 的 嗅觉 系统 在 其 生 
存 和 生殖 中 至 关 重 要 。 研 究 昆虫 的 嗅觉 系统 有 助 于 
阐明 昆虫 个 体 间 不同 昆虫 物种 间 昆虫 与 植物 间 的 
化 学 通讯 机 制 ,为 高 效 、 高 特异 性 的 昆虫 行为 调节 剂 
研发 提供 理论 基础 与 分 子 骤 标 (Plettner, 2002)。 近 
年 来 转录 组 高 通 量 测序 技术 的 快速 发 展 极 大 地 促进 
了 昆虫 功能 基因 的 研究 ( 张 棋 衣 和 豆 明 龙 , 2013 ) 。 
转录 组 测序 的 广泛 应 用 也 为 解决 昆虫 特别 是 鳞 翅 目 
昆虫 的 分 类 、 毒 理 、 发 育 、 与 寄主 互 作 以 及 非 编 码 











































































































一 些 近 缘 物种 的 嗅觉 基因 进行 进化 分 析 ,为 生物 控 
制 农林 害虫 提供 重要 的 参考 和 借鉴 。 


1 材料 与 方法 


1.1 供 试 昆 虫 

供 试 荔枝 带 蛙 虫 成 虫 采集 于 广东 省 农业 科学 院 
果树 研究 所 生态 荔枝 龙眼 园 , 在 实验 室内 通过 显 微 
镜 镜 检 区 分 荔枝 蒂 星 虫 具 雄 成 虫 ,并 盘 选 活力 充沛 的 
个 体 样 品 (雌雄 成 虫 各 5 头 ) 混合 用 于 转录 组 测序 。 
1.2 总 RNA 的 提取 与 检测 

采用 Trizol Reagent. 方法 分 别提 取 以 上 荔枝 蒂 
星 虫 混合 样品 的 总 RNA, 1% 琼脂 糖 凝 胶 电 泳 检测 
总 RNA 的 质量 。Nanodrop 分 光 光 度 计 (IMPLEN， 
CA, USA) 检 测 RNA B5 SE ( OD s0 ODs EEIE ) , 根 
据 Qubit RNA Assay Kit 说 明 对 总 RNA 浓度 进行 精 
确定 量 ; 最 后 用 Agilent 2100 ( Agilent Technologies, 
CA, USA) 精确 检测 RNA 的 完整 性 。 
1.3 文库 构建 与 测序 

取 检 测 合格 的 RNA 样品 1.5 pg, 用 带 有 Oligo 
(dT) 的 磁 珠 富 集 mRNA ,随后 将 mRNA 打 断 成 短片 
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段 。 以 mRNA 为 模板 进行 150 bp 测序 文库 的 构建 ， 
并 使 用 AMPure XP system 对 其 进行 纯化 。 最 后 将 
合成 的 cDNA 样品 在 Ilumina HiSeq'" 4000 测序 平 
台 上 进行 测序 。cDNA 文库 构建 与 测序 由 北京 诺 禾 
致 源 生物 信息 科技 有 限 公司 协助 完成 。 
1.4 数据 的 拼接 与 组 装 

测序 完成 后 获得 原始 reads ,去 除 带 接头 .ploy-N 
和 低 质 量 的 reads, 获得 clean reads。 同 时 ,计算 
Q20, Q30, GC 含量 和 重复 序列 水 平 。 采 用 Trinity 
软件 对 clean reads 进行 拼接 .过 滤 和 组 装 后 获得 得 
到 高 质量 的 unigene。 
1.5 Unigene 的 功能 注释 

分 别 以 NCBI 蛋白 数据 库 ( NCBI non-redundant 
protein sequences, Nr); NCBI 核酸 序列 数据 库 
( NCBI non-redundant nucleotide sequences, Nt ) ; £& 
白 家 族 ( protein. family, Pfam); 直系 同 源 基 因 簇 
(clusters of orthologous groups of proteins, KOG/ 
COG) ; Swiss-Prot 数据 库 (a manually annotated and 
reviewed protein sequence database ) ; 京都 基因 与 基 
因 组 百科 全 书 (KEGG ortholog database, KO) ;基因 
本 体 论 ( Gene Ontology, GO) 数据 库 做 参考 , 进行 
BLAST 程序 比 对 , 通过 gene 的 相似 性 进行 功能 
注释 。 
1.6 嗅觉 相关 基因 分 析 

根据 已 经 获得 的 功 校 带 蛙 虫 转录 组 数据 ,对 嗅 
觉 系统 相关 基因 进行 统计 分 析 。 对 NCBI. 公 共 数 据 
库 进 行 检 索 ,搜集 已 完成 全 基因 组 序列 测序 和 已 公 
开发 表 的 常见 鳞 翅 目 昆 虫 OBP 基因 的 氨基 酸 序列 ， 
将 测 得 的 荔枝 蒂 星 虫 OBP 转录 组 数据 和 己 知 的 鳞 
xH H EL: OBP 序列 进行 ClustalX 多 序列 比 对 ,对 
OBP 特有 的 6 个 保守 的 半 胱 氨 酸 残 基 模 式 进行 检 
索 分 析 。 并 使 用 MEGA7.0(Kumar et al., 2016) 的 
Neighbor-Joining 方法 ,以 bootstrap =1 000 构建 进化 
树 (Felsenstein, 1985; Saitou and Nei, 1987) 。 



































2 结果 


2.1 荔枝 蒂 星 虫 转录 组 的 测序 和 组 装 

荔枝 带 星 虫 转录 组 测序 共产 出 45 529 172 条 
reads ,去 除 低 质量 的 和 含有 接头 的 reads 以 后 ,得 到 
44 818 126 条 clean reads , 共 测 得 6.72 G MIZER 
(nucleotides, nt) 。G - C 含量 为 44. 24% , Q30 为 
92.2396 ( >90% ) (€ 1) 。 数 据 显示 测序 质量 好 ,可 


信和 度 高 。 




















表 1 Illumina HiSeq”™4000 测序 产 出 质量 统计 


Table 1 Quality of Ilumina HiSeq™ 4000 sequencing data 


原始 Reads 数 高 质量 Reads 数 ”数据 大 小 Q20 Q30 (G+C) 





Number of Number of Size of (9) (%) " 
raw reads clean reads clean bases 
45 529 172 44 818 126 6.72 G | 96.95 92.23 44.24 





Qao ,Qa0: 分 别 表示 质量 值 三 20 和 30 的 碱 基 所 占 百分比 Percentages 


of bases with the quality value =20 and 30, respectively. 
利用 Trinity 软件 对 这 些 reads 进行 组 装 , 取 每 
条 基因 中 最 长 的 转录 本 作为 Unigene ,得 到 了 68 996 
条 unigenes( 表 2) ,其 中 N50 和 平均 长 度 分 别 为 
1 444 和 762 bp。 拼 接 unigene 与 转录 本 长 度 分 布 显 
示 , 长 度 大 于 1 000 bp 的 unigene 有 14 242 条 , 占 全 部 
Unigenes 的 20. 6496 。 说 明 本 研究 中 转录 组 文库 的 测 
序 和 组 装 结 果 较 好 ,能 够 进行 后 续 生 物 信 息 学 分 析 。 
m2 荔枝 蒂 星 虫 转录 组 组 装 质量 统计 

Table 2 Assembly quality of Conopomorpha sinensis 

转录 本 Transcripts Unigenes 
数量 ”百分比 ”数量 百分比 


Number Percentage Number Percentage 











长 度 范围 (bp ) 


Sequence length 






































«301 34936 32.95 27744 40.21 

301 -500 23619 22.28 15395 22.31 

501 -1 000 21158  Á 19.95 11615 16.83 

1 001 -2 000 15188 14.33 8 214 11.91 

»2 000 11119 10.49 6 028 8.74 
总 数 Total number 106 020 68 996 
N50 长 度 N50 length (bp) 1 581 1 444 
平均 长 度 Mean length (bp) 863 762 


N50: fixi; 5096 所 有 核 苷 酸 的 最 大 Unigene KJE Length of the smallest 
contig which when added to a set of larger contigs yields at least 50% of 


the genome. 


2.2 Unigenes 的 功能 注释 

使 用 BLAST 程序 将 组 装 得 到 的 unigenes 
Ej Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KO 和 GO 
数据 库 进 行 比 对 , 获得 unigene 的 蛋白 功能 注释 信 
息 。 结 果 如 表 3 所 示 , 共 有 2 443 条 unigenes 在 以 
上 7 个 数据 库 中 同时 标注 成 功 , 占 总 unigenes 数 的 
3.54% 。 并 且 在 以 上 7 个 数据 库 中 至 少 1 条 数据 库 
注释 成 功 的 unigenes 有 22 348 条 , 占 总 unigenes 数 
的 32.39% 。 其 中 ,在 Nr 注释 成 功 的 unigenes 的 数 
量 最 多 (27.01% ) ,其 后 依次 是 G0(22. 58% ) , Pfam 
(22. 4196 ) , Swiss-Prot ( 18. 88% ) , KOG ( 14. 33% ) ， 
KEGG(10. 53% ) 和 Nt (7. 4396) , Unigene 注释 到 
Nr 数据 库 中 的 物种 分 布 图 显示 (图 1) : DRR 
Plutella xylostella 所 占 比 例 最 高 (34. 196 ) ,其 后 依次 
是 家 看 Bombyx mori(25.9% ) ,帝王 蝶 Danaus plexippus 
(17. 996) , THE iR UIS Papilio xuthus (1. 496 ) , i «3 Wr 
Acyrthosiphon pisum(1. 496 ) ,其 他 物种 占 19.3% 。 
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A3 荡 枝 蒂 星 虫 基因 注释 成 功率 统计 
Table 3 Successful annotation rate of genes of 


Conopomorpha sinensis 











注释 数据 库 独立 基因 注释 百分比 
Annotated database Unigenes Percentage 
Nr 18 637 27.01 
Nt 5 130 7.43 
KEGG 7 272 10.53 
Swiss-Prot 13 028 18.88 
Pfam 15 467 22.41 
GO 15 585 22.58 
KOG/COG 9 889 14.33 
VE sols MERI 
Mecum al dabas 0249 S 
14S crx 
md u ps ee d 
总 的 unigenes Total unigenes 68 996 100 
物种 分 布 Species distribution 
34.196 






25.9% 


19.3% 


17.9% 1.4% 


1.4% 
小 菜 蛾 Plutella xylostella 
dd Bombyx mori 

I EHE Danaus plexippus 
WEY Acyrthosiphon pisum 
HIRR Papilio xuthus 


其 他 Others 


Unigenes 在 Nr 数据 库 中 的 物种 分 布 图 


Species distribution of unigenes in Nr database 








图 1 
Fig. 1 








2.3 Unigenes 的 功能 分 类 

基于 GO 数据 库 对 功 校 带 蝗虫 的 unigenes 进行 
功能 分 类 统计 ,可 以 从 宏观 上 认识 荔枝 蒂 星 虫 的 基 
因 功 能 分 布 特征 。 对 荔枝 带 旺 虫 转录 组 unigene yt 
行 GO 分 析 发 现 ,在 GO 注释 的 15 585 条 unigene 
中 ,注释 到 生物 学 过 程 (biological process) 的 基因 最 
多 (195 466) ,其 次 是 细胞 组 分 (cellular component) 
(58 600), 分 子 功 能 (molecular function ) 的 最 少 
(35 532) 。 在 47 个 亚 类 中 ,以 细胞 进程 和 代谢 过 程 
在 生物 过 程 中 的 比例 较 高 ,细胞 和 细胞 成 分 在 细胞 
组 分 中 所 占 比例 较 高 ,蛋白 结合 和 催化 活性 在 分 子 
功能 分 类 中 占有 较 高 比例 (图 2)。 

结合 KEGG Zt EXT zx scs sk HT unigenes 可 
能 参与 或 涉及 的 代谢 途径 进行 分 析 , 结果 表明 ， 
7 272 条 unigenes 参与 了 5 类 代谢 通路 分 支 , 即 细胞 
过 程 ( cellular processes ), 环境 信息 处 理 


(environmental information processing) , 遗传 信息 处 

















理 ( genetic information processing ) ,代谢 (metabolism ) 
和 有 机 系统 (organismal systems ) ,267 个 代谢 通路 。 
其 中 注释 最 多 的 6 个 代谢 通路 依次 为 :核糖 体 
(302) ,内 质 网 蛋白 加 工 (242 ) ,剪接 体 (193 ) , 碳 代 
谢 (188 ) ,氧化 磷酸 化 (182) , 味 吟 代谢 (168)( 表 4) 。 
2.4 嗅觉 相关 基因 分 析 

通过 对 7 大 公共 数据 库 进 行 同 源 比 对 和 基因 注 
释 ,100 个 蔓 术 蒂 星 虫 嗅觉 相关 基因 在 功 枝 人 蒂 星 虫 
转录 组 中 获得 注释 。 主 要 包括 : 41 个 气味 结合 蛋 
白 (odorant binding proteins, OBPs) 基因 ,11 个 化 学 
感受 蛋白 (chemosensory proteins, CSPs) 基因 ,45 个 嗅 
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Cellular Component 


图 2 G0 分 类 图 


Molecular Function 





Fig. 2 Go categorization 
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表 4 KEGG 注释 上 unigenes 最 多 的 6 个 代谢 通路 
Table 4 Six unigene pathway hierarchies with the most 
annonated unigenes of Conopomorpha sinensis on KEGG 


代谢 通路 基因 数 








e Pathway hierarchy Number of genes 
ko03010 核糖 体 Ribosome 302 
ft uil ade in endoplasmic reticulum 242 
ko03040 剪接 体 Spliceosome 193 
ko01200 碳 代 谢 Carbon metabolism 188 
ko00190 氧化 磷酸 化 Oxidative phosphorylation 182 
ko00230 IE [Cif Purine metabolism 168 


觉 受 体 (olfactory receptors, Ors) 基因 和 3 个 离子 型 
ZAA (ionotropic receptors, IRs) 基因 ( 表 5)。 其 中 ， 
有 11 个 气味 结合 蛋白 基因 ,11 个 化 学 感受 蛋白 基 
因 和 3 个 离子 型 受 体 基因 在 Nr 数据 库 中 获得 注释 。 
m5 荔枝 蒂 蛙 虫 转录 组 中 与 嗅觉 相关 基因 统计 


Table 5 Statistics of unique transcripts associated with 





olfactory of Conopomorpha sinensis 


被 注释 基因 数 Nr 注释 数 
















































































嗅觉 相关 基因 
人 HX Number of annotated Nr annotation 

Olfaction-related genes : 

unigenes number 

气味 结合 蛋白 基因 新 f 
Odorant binding protein ( OBP) genes 
化 学 感受 蛋白 基因 
Chemosensory protein ( CSP) genes 
嗅觉 受 体 基因 
Olfactory receptor (Or) genes 
离子 型 受 体 基 ， ， 


Ionotropic receptor (IR) genes 


通过 Novofinder 搜索 unigene 数据 库 中 荔枝 蒂 
旺 虫 气味 结合 蛋白 OBPs 基因 ,我 们 共 获 得 39 条 
unigene 编码 功 校 带 蝗 虫 气味 结合 集 白 基因 序列 。 
其 中 包含 37 2& 2x osi sk d, OBP 基因 ,命名 为 “C. 
sinensis OBPx" ;2 条 荔枝 蒂 星 虫 GOBP 基因 ,命名 为 
"C. sinensis GOBPx" , 3X8 it 5j E, An XH EL EG E 
OBPs 的 蛋白 序列 比 对 分 析 ,我们 鉴定 出 12 个 基因 
的 蛋白 序列 拥有 OBP 典型 的 6 个 保守 半 胱 氨 酸 模 
式 。 系 统 进化 分 析 发 现 (图 3) , 23 Bo kk t COBP 
分 文 聚 类 较为 集中 ;大 部 分 C. sinensis OBP 都 至 少 与 
一 个 鳞 翅 目 OBP 同 源 物 聚 类 在 一 支 , 且 这 些 序 列 与 
NCBI 上 登录 的 序列 拥有 75% 以 上 的 氨基 酸 序列 一 
致 性 。 对 于 鳞 翅 目 不 同 种 的 昆虫 ,也 出 现 一 定 的 上 聚 
类 关系 ,如 荔 校 带 旺 虫 中 预测 的 OBPs (C. sinensis 
OBPI2 和 OBP19; C. sinensis OBP21, OBP23, 
OBP28 和 OBP37) ,/ 3E; OBPs( P. xylostella OBP3 
和 P. xylostella PBP1) PIE X$ OBPs(O. nubilalis 



































PBP3) 分 别 在 进化 树 中 形成 一 个 独立 的 小 分 文 , 表 
明 昆 虫 种 类 分 化 后 也 出 现 了 各 自 OBPs 的 分 化 。 





3 讨论 





随 着 高 通 量 测序 技术 的 发 展 ,转录 组 分 析 逐 渐 
成 为 非 模式 昆虫 中 发 掘 功能 基因 的 一 种 有 效 手 段 
(Zou et al., 2008; Vogel et al., 2011) ,并 趋向 于 解 
决 有 害 生 物 防治 `. 近 缘 种 遗传 关系 与 进化 等 更 加 复 
杂 的 实际 问题 (Yang et al., 2014) 。 本 研究 首次 利 
用 RNA-Seq 技术 对 荔枝 蒂 星 虫 的 转录 组 进行 了 测 
序 和 功能 分 析 。 研 究 共 获 得 68 996 条 质量 较 高 的 
unigenes。 通 过 数据 库 分 析 和 基因 功能 注释 ,深度 发 
掘 并 获得 100 4 5j zy dotis Bk rb obse 4H OS B SE bs 
基因 。 

基因 功能 注释 结果 显示 ,通过 Blast, [B] UR TETE 
索 获 得 22 348 条 unigenes (32. 3996 ) 成 功 注释 , 仍 有 
46 648 条 unigenes (60. 61% ) 未 被 注释 。 这 可 能 是 
unigene 的 长 度 直接 影响 着 基因 功能 注释 的 结 
unigene 较 短 而 未 与 公共 数据 库 中 的 序列 比 对 上 。 
也 可 能 是 非 编码 序列 或 者 是 新 的 基因 (Konishi and 
Ohnishi, 2006; Hou et al., 2011), 。 此 外 ,基因 数据 
库 中 的 参考 序列 的 数量 是 决定 转录 组 数据 注释 的 另 
一 大 因素 。 例 如 , 理 豆 蚜 的 基因 组 为 半 翅 目 昆 虫 的 
序列 注释 提供 了 极 大 便利 ( Eisen, 2010) ;已 有 的 烟 
AL Bemisia tabaci 和 褐飞虱 Nilaparvata lugens 的 基 
因 序列 信息 也 为 其 他 粉 天 和 飞 剧 类 昆虫 的 分 子 生物 
学 研究 提供 了 便利 (Wang et al., 2010; Xue et al., 
2014) ;家 乔 作 为 鳞 运 目 昆虫 研究 的 模式 昆虫 ,其 基 
因 组 全 序列 的 测定 也 为 鳞 翅 目 昆 虫 的 比较 生物 学 和 
功能 基因 组 研究 奠定 了 基础 (Xia et al., 2004) , H 
前 ,高 通 量 测序 已 在 帝王 蝶 ( Zhan et al., 2011) 、 小 
菜 蛾 (You et al., 2013) AEI Cnaphalocrocis 
medinalis (Li et al., 2013) , 4 ZX T IR Spodoptera 
litura( Gu et al., 2013) 等 多 种 鳞 翅 目 非 模式 昆虫 中 
开展 ( 张 传 溪 , 2015 ) 。 因 此 ,对 荔枝 带 蛙 虫 的 高 通 
量 测 序 分 析 将 推动 鳞 翅 目 非 模式 昆虫 研究 的 全 面 
发 展 。 

利用 GO 数据 库 对 范 校 带 旺 虫 unigenes 进行 基 
因 功 能 分 类 ,可 以 全 面 描 述 生 物体 中 基因 和 基因 产 
物 的 属性 。KEGG 数据 库 分 析 表 明 , 功 枝 带 旺 忠 
7 272 条 unigenes 定位 的 267 个 代谢 途径 分 支 分 别 
参与 到 荔枝 蒂 星 虫 体内 的 核糖 体 代 谢 、 碳 水 化 合 物 
代谢 、 氧 化 磷酸 化 等 过 程 中 。 其 中 ,包含 unigenes 最 
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图 3 ”基于 氨基 酸 序列 构建 的 荔枝 蒂 性 虫 嗅觉 受 体 CsOBP 与 其 鳞 翅 目 近 缘 物种 OBP 的 系统 发 生 关系 

Fig. 3 Phylogenetic tree of CsOBPs with OBPs from the known closely related species of Lepidoptera based on amino acid sequences 

OBPs 的 GenBank 登录 号 The GenBank accession numbers of OBPs; B. mori GOBPI ; NM, 001044031. 1; B. mori GOBP2; NM. 001044033. 1; B. 

mori BOP1; AB353892.1; B. mori BOP2; NM. 001163900. 1; B. mori BOP3; NM 001160191.1; B. mori BOP4: AB508819. 1; B. mori BOPS : 
AB508818.1; B. mori BOP6; AB353894.1; S. litura GOBP1; KP331514. 1; S. litura GOBP2; KP331515. 1; S. litura OBP1: KP331516. 1; S. 

litura OBP2; KP331517. 1; S. litura OBP3; KP331518. 1; S. litura OBP4: KP331519. 1; S. litura OBPS; KP331520. 1; S. litura OBP6: 
KP331521.1; S. litura OBP7; KP331522. 1; S. litura OBP8; KP331523. 1; S. litura putative OBP16: KT261662. 1; S. litura putative OBP21: 
KT261667.1; S. litura putative OBP9; KT261655. 1; S. litura putative OBP24; KT261670. 1; P. xylostella GOBP2; EU368115. 2; P. xylostella 
GOBPI ; EU368114.1; P. xylostella OBP1; JNO88193. 1; P. xylostella OBP2; JNO88194. 1; P. xylostella OBP3; JNO88195. 1; P. xylostella OBP4 ; 
JN088196.1; P. xylostella OBP5; JNO88197.1; P. xylostella PBP1; FJ201994. 1; P. xylostella PBP2; JX308238.1; O. nubilalis PBP1; GU828019. 

1; O. nubilalis PBP2; GU828020. 1; O. nubilalis PBP3: GU828021.1; O. nubilalis PBP4; GU828022. 1; O. nubilalis PBP5; GU828023.1; S. 

exigua GOBP1; KF840411.1; S. exigua OBP1; JX962784.1; S. exigua OBP2; JX962785.1; S. exigua OBP3; JX962786. 1; S. exigua OBP4: 
JX962787.1; S. exigua OBPS; JX962788. 1; S. exigua PBPl; EF052846. 1; S. exigua PBP2; EF052849. 1; P. xuthus GOBPI-like; XM _ 
013307375. 1; P. xuthus GOBP2-like; XM. 013307373.1; D. kikuchii GOBP1; KC783391.1; D. kikuchii GOBP2; KC783387.1; D. kikuchii OBP1 ; 
KF487589.1; D. kikuchii OBP2; KF487590. 1; A. segetum GOBP1: DQ838577.1; A. segetum GOBP2; DQ838579.1; C. suppressalis GOBPI ; 
EU825768.1; C. suppressalis GOBP2: DQ447635. 1; C. suppressalis OBP1; JNO88201. 1; M. sexta GOBP2; AF323972. 1; M. sexta PBPl: 
AF323972.1; A. ipsilon GOBP1; DQ838578.1; A. ipsilon GOBP2; JQ822244. 1; A. ipsilon OBP1; JX863689. 1. 以 bootstrap = 1 000 构建 进化 树 。 


Bootstrap support values presented as percentages were based on 1 000 replicates. 





























多 的 是 核糖 体 代 谢 通路 (ko03010 ) , 共 302 条 ,这 可 合成 部 分 核糖 体 蛋 白质 ,并 形成 各 种 酶 (Guertin et 
能 是 生物 体 的 生长 发 育 和 行为 调控 与 核糖 体 数 量 和 al., 2006) ,然后 这 些 酶 进一步 促使 新 的 mRNA 生 
代谢 相关 导致 的 结果 ,核糖 体 蛋白 基因 作为 核糖 体 成 ,合成 更 多 的 核糖 体 蛋 白质 ,从 而 导致 核糖 体形 
的 重要 组 成 部 分 ,对 和 蛋白质 的 合成 举足轻重 ,从 而 影 ”成 。 本 研究 为 进一步 大 量 挖掘 荔 杖 带 星 虫 生理 调控 
响 生 物 的 整个 生命 活动 (Uechi et al., 2001), H ” 中 的 重要 功能 基因 ,开展 荔枝 带 蛙 虫 靶 标 基 因 克 隆 
带 蛙 虫 在 生命 活动 过 程 中 ,体内 的 mRNA 可 能 指导 ”及 功能 验证 等 提供 了 基础 数据 。 
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嗅觉 系统 对 于 昆虫 的 生命 活动 至 关 重 要 ( Field 
et al., 2000) 。 在 长 期 进化 过 程 中 ,昆虫 形成 了 高 度 
灵敏 的 嗅觉 系统 ,并 以 此 来 感受 自然 环境 中 的 各 种 
化 学 信息 ,产生 相应 的 生理 和 行为 反应 ,如 寻找 配 
偶 、 食 物 源 和 栖息 场所 、 产 卵 行为 .逃避 危险 或 者 不 
适合 的 栖息 地 和 寄主 等 (Karban and Baldwin, 1997; 
Will and van Bel, 2006; Walling, 2008), 。 因 此 , 对 
昆虫 的 嗅觉 系统 进行 研究 继而 对 嗅觉 发 生 过 程 进行 
干扰 是 防治 害虫 的 一 条 有 效 途 径 (Plettmner,，2002 ) 。 
本 研究 首次 从 荔 校 带 旺 虫 转录 组 数据 库 中 发 据 出 
100 个 嗅觉 相关 基因 ,包含 :45 个 嗅觉 受 体 CsORs 
基因 、3 个 离子 型 受 体 CsIRs 基因 、41 个 气味 结合 
EA CsOBPs 基因 和 11 个 化 学 感受 蛋白 CsCSPs 基 
。 前 人 对 其 他 鳞 翅 目 昆 虫 基因 组 和 转录 组 测序 分 
Jr ,也 发 现 了 不 同 种 类 和 数量 的 嗅觉 相关 基因 。 如 ; 
小 菜 蛾 成 虫 触角 转录 组 中 ,鉴定 获得 110 个 嗅觉 相 
关 基 因 ,包含 54 个 嗅觉 受 体 基 因 、16 个 离子 型 受 体 
基因 25 个 气味 结合 蛋白 基因 和 15 个 化 学 感受 和 蛋 
白 基 因 i — fb Chilo suppressalis 成 里 触角 转录 组 
的 分 析 共 识别 出 47 个 CsupOR 基因 、20 个 CsupIR 
基因 、26 个 CsupOBP 基因 和 21 个 CsupCSP 基因 
(Cao et al., 2014) 。 相 比 与 其 他 鳞 却 目 转录 组 研 
究 ,我 们 获得 荔枝 带 蝗虫 嗅觉 基因 数量 处 于 合理 范 
畴 ,但 在 离子 型 受 体 基因 和 气味 结合 蛋白 基因 也 存 
在 数量 的 差异 性 。 研 究 进 一 步 与 家 天、 斜纹 夜 蛾 、 小 
菜 蛾 、 玉 米 晨 等 鳞 怒 目 近 缘 物 种 气味 结合 蛋白 基因 
联合 分 析 发 现 , 与 荔枝 蒂 蛙 虫 直 系 同 源 的 气味 结合 
蛋白 基因 18 组 ;同时 也 存在 一 批 仅 在 范 术 带 旺 虫 中 
存在 的 气味 结合 蛋白 基因 。 这 些 差 异 和 独 有 的 基因 
可 能 与 荔枝 带 旺 虫 物种 生境 中 特有 的 化 学 物质 相 
关 , 并 与 它 隐 蔽 性 强 和 寄主 专 一 特性 息息相关 。 
此 深入 研究 荔枝 蒂 星 虫 嗅觉 识别 的 分 子 机 制 , 将 为 
高 效 .高 特异 性 的 昆虫 行为 调节 剂 研发 提供 理论 基 
础 与 靶 标 。 此 外 ,通过 未 梳 蒂 星 虫 转录 组 分 析 , 更 多 
具备 特异 性 高 效 性 安全 性 的 昆虫 防治 靶 标 将 被 发 
现 ,为 无 公害 荔 校 生产 提供 安全 、 低 残留 的 害虫 控制 
措施 。 
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